Compromise and Synchronization 
in Opinion Dynamics 
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We discuss two models of opinion dynamics. First we present a brief review of the Hegselmann 
and Krause (HK) compromise model in two dimensions, showing that it is possible to simulate the 
dynamics in the limit of an infinite number of agents by solving numerically a rate equation for a 
continuum distribution of opinions. Then, we discuss the Opinion Changing Rate (OCR) model, 
which allows to study under which conditions a group of agents with a different natural tendency 
(rate) to change opinion can find the agreement. In the context of the this model, consensus is 
viewed as a synchronization process. 
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Since the behavioral revolution and the birth of cy- 
bernetics, the so called 'soft' social sciences have emu- 
lated both the intellectual and methodological paradigms 
of the 'strong' natural sciences [[J. The certainty and 
stability of the Newtonian paradigm has represented for 
decades the cornerstone of sciences like psychology, econ- 
omy and sociology, which have been largely inspired 
by classical mechanics and statistical thermodynamics. 
Clearly this trend has continued since quantum mechan- 
ics, chaos and complexity revolutions have leaded to a re- 
consideration of the relevance of the Newtonian paradigm 
to all natural phenomena. 

In the last years, new disciplines such as econophysics 
and sociophysics have largely demonstrated the power of 
agent-based computational models in simulating complex 
adaptive systems (financial markets, cultural evolution, 
social structures, voter communities) in which large num- 
bers of individuals are involved in massively parallel local 
interactions 0,0 ■ In agent-based models, individuals are 
modelled as autonomous interacting agents with a vari- 
able degree of internal complexity and numerical simula- 
tions represent computational experiments to study the 
evolution of a given social system under controlled con- 
ditions. Of course in many cases the individual cogni- 
tive behavior is oversimplified, as for example in opinion 
dynamics models where human opinions are reduced to 
integer or real numbers [JjQ • In more complicated mod- 
els, individuals are simulated by means of simple neural 
networks or associative memories. Although in this case, 
many of the simplifications adopted are somehow unreal- 
istic. On the other hand, also the Kepler's laws assump- 
tion of Earth as a point-mass was not realistic at all, but 
for the purpose of describing celestial motion it turned 
out very successful. Furthermore, the aim of agent-based 
simulations is to provide information on averages over 
many people, and not on the fate of a specific person. In 
this sense, despite of their simplicity, these models seem 
to work very well. For example, the Sznajd model pre- 
diction of the distribution of votes among candidates in 
Brazilian and Indian elections is encouraging, although 
the model is not able to predict the number of votes one 



specific candidate gets in one specific election 0- 
In the first part of this paper we discuss one of the well 
known models of opinion dynamics, the so called com- 
promise model of Hegselmann and Krause Q, showing 
some recent results about its continuum version in a two- 
dimensional opinion space 0. In the second part of the 
paper we discuss a new perspective in opinion dynamics 
based on the suggestion of a possible role of synchroniza- 
tion in opinion formation. By means of the so called 
Opinion Changing Rate model, a modified version of the 
Kuramoto model adapted to a social context, we study 
under which conditions a group of agents with a differ- 
ent natural tendency (rate) to change opinion can find 
agreement 8] . 



I. DISCRETE AND CONTINUUM OPINION 
DYNAMICS IN THE 2- VECTOR HK 
CONSENSUS MODEL 

For the sociologist Robert Axelrod "culture" is mod- 
elled as an array of features, each feature being specified 
by "traits" , which are expressed by numbers. The num- 
ber of features or dimensions is nothing but the number 
of components of a vector, and two persons interact if 
and only if they share at least one common feature (i.e. 
the same value of the corresponding vector component) 
H • In this model, two persons are culturally closer the 
more features they have in common, and the number of 
these common features is, in turn, related to the proba- 
bility for the two individuals to interact. Starting from 
the Axelrod model, several simple agent-based models of 
opinion formation have been devised, mostly by physi- 
cists HSU. 

In general, a typical scalar opinion formation model 
starts by assigning randomly a real number (chosen in 
a given interval) to every agent of the system. Then the 
dynamics starts to act, and the agents rearrange their 
opinion variables, due to their interactions. At some 
stage, possibly, the system reaches a configuration which 
is stable under the dynamics. This final configuration 
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may represent consensus, when all agents share the same 
opinion, polarization, when there are two main clusters 
of opinions ("parties"), or fragmentation, when several 
opinion clusters survive. However, a discussion between 
two persons is not simply stimulated by their common 
view/preference about a specific issue, but it in general 
depends on the global affinity of the two persons, which 
is influenced by several factors. So, for a more realistic 
modelling of opinion dynamics, one should represent the 
opinions/attitudes like vectors (as in the Axelrod model), 
and not like scalars. In this section we will focus on the 2- 
vector version of the Hegselmann and Krause (HK) com- 
promise model and we will show that it is possible to 
simulate the discrete opinion dynamics in the limit of an 
infinite number of agents by solving numerically a rate 
equation for a continuum distribution of opinions Q • 
The HK model is based on the concept of bounded 
confidence, i.e. on the presence of a parameter e, 
called confidence bound, which expresses the compatibil- 
ity among the agents in the opinion space. If the opinions 
of two agents i and j differ by less than e, their positions 
are close enough to allow for a discussion, which even- 
tually leads to a change in their opinions, otherwise the 
two agents do not interact with each other. The physical 
space occupied by the agents living in a society or a com- 
munity can be modelled as a graph, where the vertices 
represent the agents and the edges relationships between 
agents. So we say that two agents can eventually talk 
to each other if there is an edge joining the two corre- 
sponding vertices (in graph language, if the two vertices 
are neighbors). In the following we will consider only 
the general case of a society where everybody talks to 
everybody. The dynamics of the HK model is usually 
simulated by means of Monte Carlo (MC) algorithms. 
One chooses at random one of the agents and checks how 
many of its neighbors (in the physical space) are com- 
patible, i.e. lie inside the confidence range in the opinion 
space. Next, the agent takes the average opinion of its 
compatible neighbors. The procedure is repeated by se- 
lecting at random another agent and so on. The type 
of final configuration reached by the system depends on 
the value of the confidence bound e. For a scalar opinion 
space [0, 1] it has been shown that consensus is reached 
for e > e c , where the critical threshold e c is strictly re- 
lated to the type of graph adopted to model society: ac- 
tually, it can take only one of two possible values, e c ~ 0.2 
and 0.5, depending on whether the average degree of the 
graph (i.e. the average number of neighbors) diverges, 
as in our case of a completely connected graph, or stays 
finite when the number of vertices goes to infinity |1C)| . 
On the other hand, the 2-vector HK model on a com- 
pletely connected graph is much less studied than the 
1-dimensional version. In this case the opinion space is 
represented by the points (x, y) of a bidimensional mani- 
fold, that in general is a square [0, 1] x [0, 1] and the con- 
fidence range is a circle whose radius is the confidence 
bound e. 

In Fig.l we plot a sequence of snapshots of the opinion 
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FIG. 1: Sequence of snapshots of the 2D squared opinion 
space for the 2-vector discrete HK model with N = 2000 
agents and e = 0.25. The points in the upper-left panel rep- 
resent different randomly distributed opinions at t=0. In the 
other panels, where we show successive (but not consecutive) 
time steps of a Monte Carlo simulation with simultaneous 
update, opinions merge together in different clusters. Finally, 
in the lower-left panel (t=12) the consensus is reached and 
all the opinions occupy the same position. In the second-last 
panel, the number of opinions concentrated in the bottom 
clusters is also indicated (see text). 




FIG. 2: Sequence of snapshots of the 2D circular opinion space 
for the 2-vector discrete HK model with N = 2000 agents 
and e = 0.25. The dynamics is the same than in Fig.l. The 
different shape of the space influences the dynamical evolution 
of the opinions and consensus is reached in a shorter time 
(t=6) than in the case of Fig.l. 



space for a discrete configuration of N = 2000 agents 
and a value e = 0.25. Each point in the first snapshot 
(upper-left panel) represents the opinion of one agent at 
time t=0. Then the system evolves by means of simulta- 
neous updates (i.e. all the opinions are updated at each 
MC time step) merging the opinions in bigger and big- 
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ger clusters until a stationary state is reached. After 12 
time steps (lower-right panel) consensus is fully obtained 
and all the opinions lie on the same cluster. More in 
general, from extended numerical simulations it results 
yfl that the consensus threshold for the discrete 2-vector 
HK model is e c ~ 0.24, a value slightly greater than that 
one found for the scalar model with the same topology. 
This value tends to the value e c ~ 0.23 when the number 
of agents grows. 

If we look at the basis of the triangle in the second-last 
snapshot of Fig.l, just before reaching final consensus, we 
can see that the two big clusters at the vertices, made by 
around 500 opinions and at a reciprocal distance greater 
than the confidence bound (thus a-priori not interact- 
ing), are in contact only by means of a small cluster of 21 
agents. Such a phenomenon is very frequent in the Monte 
Carlo simulations of the HK model when the system is 
approaching consensus (in both one and two dimension) ; 
indeed, almost always consensus is reached only because 
of this phenomenon. This models an important feature 
of real social networks, i.e. the existence of the so called 
connectors which play the role of a bridge between oth- 
erwise not interacting social grou ps, thus ensuring the 
cohesiveness of the entire network |ll| . 
Another peculiar feature of the HK model, clearly visi- 
ble in the upper snapshots of Fig.l, is the fact that the 
dynamics always starts to act from the edges of the opin- 
ion space, where the opinion distribution is necessarily 
inhomogeneous, so that it is essentially the shape of the 
opinion space which rules the symmetry of the resulting 
cluster distribution. In order to better appreciate this ef- 
fect, we plot in Fig. 2 the temporal evolution of the same 
system of Fig.l, but with a circular opinion space. In 
this case, even if the final configuration is the same as 
before, the resulting dynamics is different and, for the 
same value of the confidence bound, consensus is reached 
more quickly (six MC time steps with simultaneous up- 
date) due to the greater symmetry of the opinion space. 
The circular symmetry has a remarkable effect also on 
the consensus threshold, that in this case tends to that 
of the corresponding scalar HK model, i.e. e ~ 0.2. In 
a recent paper it has been shown that the 2-vector 
HK model on a completely connected graph and with a 
squared opinion space can be described by means of a 
rate equation for a continuum distribution of opinions 
P(x,y,t). The rate equation can be solved numerically, 
as already done for the scalar compromise model of Def- 
fuant et al. ^3)- The advantages of the evolution equa- 
tion over discrete Monte Carlo simulations are that one 
can directly deal with a system with an arbitrarily large 
number of agents, and that the final cluster configura- 
tions for a continuum distribution are much more sym- 
metric and regular, thus allowing a better resolution of 
the progressive merging of opinion clusters. 
This is clearly visible in Fig. 3, where we show the final 
configurations of the opinion space for several values of 
the confidence bound. In this case the squared (x, y) 
opinion space has been reduced to a grid of 100 x 100 























I. 

" 


« ft « e • o • 

ft ft ft ft o • 

a » o o a • • 
o o « « ft « « 

ft ft O O ft ft ft 
ft ft ft ft ft ft 
ft ft ft ft ft ft ft 




i. 


• • • * 

• • • • 

• • • • 

• • • ♦ 




i. 
■ 


ft O ft 

ft ft © 
ft ft ft 
































o e 
■ 

e e 




i. 


• * 

• * 




L 


o 











FIG. 3: Final configurations of the 2-vector HK model with a 
continuum distribution of opinions. From top left to bottom 
right: 6=0.08,0.10,0.15; 0.20,0.22,0.24. 



bins. All the simulations start from a flat distribution 
P(x,y,t — 0) = const and the dynamics runs until the 
distribution P(x 1 y,t) reaches a stationary state for a 
given value of the confidence bound 0. 
As one can see that, for small value of e, a regular lattice 
of clusters appears, with a squared shape inherited by the 
shape of the opinion space (as happened for the discrete 
HK model). Going on, for greater values of e, one can 
observe the progressive merging of the pairs of clusters 
with reciprocal distance less than the confidence bound 
radius. Finally, above the critical threshold e c ~ 0.23, 
consensus is completely reached. This result confirms 
the threshold value found with the MC simulations for 
a discrete dynamics of opinions in the limit of a large 
number of agents and encourages further applications of 
the rate equation technique to other opinion formation 
models. 



II. THE OPINION CHANGING RATE MODEL: 
A ROLE FOR SYNCHRONIZATION IN OPINION 
FORMATION 

Most of the opinion formation models, as for example 
the HK model presented in the previous section, have the 
limitation of not taking into account the individual incli- 
nation to change, a peculiar feature of any social system. 
In fact, each one of us changes ideas, habits, style of life 
or way of thinking in a different way, with a different ve- 
locity. There are conservative people that strongly tend 
to maintain their opinion or their style of life against 
everything and everyone. There are more flexible peo- 
ple that change ideas very easily and follow the current 
fashions and trends. Finally, there are those who run 
faster than the rest of the world anticipating the others. 
These different tendencies can be interpreted as a contin- 
uous spectrum of different degrees of natural inclination 



20 JO 
time 



■a 



10 

time 



61) 



FIG. 4: OCR model for N = 1000 and K = 1.0 (incoherent 
phase). In panels (a) and (b) we report, for each time step, 
1000 points corresponding respectively to the opinions and 
opinion changing rates of the N agents, while in panels (c) and 
(d) we show, respectively, the final distribution of opinions 
and the time evolution of the order parameter R. 



to changes. 

In a recent paper Q we have showed how such an in- 
dividual inclination to change, differently distributed in 
a group of people, can affect the opinion dynamics of 
the group itself. If we switch from the question: "Could 
agents with initial different opinions reach a final agree- 
ment ?" into the more realistic one " Could agents with 
a different natural tendency to change opinion reach a 
final agreement ? " , we can introduce a new concept, the 
natural opinion changing rate, that is very similar to the 
characteristic frequency of an oscillator. In such a way, 
we can treat consensus as a peculiar kind of synchroniza- 
tion (frequency locking), a phenomenon which has been 
very well studied in different contexts by means of the 
Kuramoto model fl3|. 

The Kuramoto model of coupled oscillators is one of the 
simplest and most successful models for synchronization. 
It is simple enough to be analytically solvable, still re- 
taining the basic principles to produce a rich variety of 
dynamical regimes and synchronization patterns. The 
dynamics of the model is given by 

K N 

6 i (t)=u H + — Y,sin(0 j -8 i ) i = l,...,N (1) 

i=i 

where 6i(t) is the phase (angle) of the ith oscillator at 
time t (— 7r < 9i(t) < 7r), while u>i is its intrinsic frequency 
randomly drawn from a symmetric, unimodal distribu- 
tion g{to) with a first moment loq (typically a Gaussian 
distribution or a uniform one) . These natural frequencies 
uji are time-independent. The sum in the above equation 
is running over all the oscillators so that this is an exam- 
ple of a globally coupled system. The parameter K > 
measures the coupling strength in the global coupling 
term. For small values of K, each oscillator tends to run 



FIG. 5: OCR model for N = 1000 and K = 2.0 (partially 
coherent phase). Same quantities as in Fig. 4. 



independently with its own frequency, while for large val- 
ues of K, the coupling tends to synchronize (in phase and 
frequency) the oscillator with all the others. Kuramoto 
showed that the model, despite the difference in the nat- 
ural frequencies of the oscillators, exhibits a spontaneous 
transition from incoherence to collective synchronization, 
as the coupling strength is increased beyond a certain 
threshold K c [lj . 

The existence of a critical threshold for synchroniza- 
tion in Kuramoto model is very similar to the consensus 
threshold found in the majority of the opinion formation 
models. Of course, at variance with the phases in the 
Kuramoto model, in a model for opinion dynamics we do 
not need periodic opinions nor limited ones: in fact, the 
opinions have a very general meaning and can represent 
the style of life, the way of thinking or of dressing etc. 
Thus we do not consider periodic boundary conditions 
and we assume Xi g] — oo, +oo[ Vi = 1, N . The dy- 
namics of the OCR model is governed by the following 
set of differential equations 

K N 

Xi(t) = u>i+—^2asin(xj-Xi)e a ^ 3 2=1 i = l,...,N 
j=i 

(2) 

Here Xi(t) is the opinion (an unlimited real number) 
of the ith individual at time t, while represents the 
so called natural opinion changing rate, i.e. the intrin- 
sic inclination, or natural tendency, of each individual 
to change his opinion (corresponding to the natural fre- 
quency of each oscillator in the Kuramoto model). As 
in the Kuramoto model, also in the OCR model the w's 
are randomly drawn from a given symmetric, unimodal 
distribution g(u>) with a first moment luq. Usually a uni- 
form distribution centered at ujq is used. In this way 
one simulates the fact that in a population there are: 1) 
conservative individuals, that naturally tend to change 
their opinion very slowly, and thus are characterized by 
a value of u>i smaller than loq\ 2) more flexible people, 
with LJi ~ luq, that change idea more easily and follow 
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FIG. 7: OCR model for N — 100 and a value of K increasing 
with a constant rate from 0.1 up to 10.1. Same quantities as 
in the previous figures. 



the new fashions and trends; 3) individuals with a value 
of LOi higher than uq, that run faster than the others in 
suggesting new ideas and insights. 

In the equation J5J) K, as usual, is the coupling strength. 
The exponential factor in the coupling term ensures that, 
for reciprocal distance higher than a certain threshold, 
tuned by the parameter a, opinions will no more influ- 
ence each other. Such a requirement is inspired by the 
confidence bound concept discussed in the previous sec- 
tion. (Please note that, due to a misprinting, there is an 
a factor missing in the coupling term of Eq.(7) of ref. 8]). 
At this point we can study the opinion dynamics of the 
OCR model by solving numerically the set of ordinary 
differential equations J2J) for a given distribution of the 
ui's (natural opinion changing rates) and for a given cou- 
pling strength K . In particular, we want to find out if, as 
a function of K , there is a transition from an incoherent 
phase, in which people change opinion each one with his 
natural rate Wj, to a synchronized one in which all the 
people change opinion with the same rate and share a 
common social trend, a sort of 'public opinion'. In order 
to measure the degree of synchronization of the system 
we decided to adopt an order parameter R{t) related to 
the standard deviation of the opinion changing rate Xj (t) 

and defined as R(t) = 1 - X^Life (t) - X(t)) 2 , 

where X(t) is the average over all individuals of Xj(t). 
It is easy to see that R = 1 in the fully synchronized 
phase, where all the agents have exactly the same opinion 
changing rate (and very similar opinions) , while R < 1 in 
the incoherent or partially synchronized phase, in which 
the agents have different opinion changing rates and dif- 
ferent opinions. The numerical simulations have been 
performed typically with N=1000 agents and with an 
uniform distribution of the initial individual opinions 
Xi(t = 0) in the range [-1,1]. The natural opinion chang- 
ing rates Wj are taken from a uniform distribution in the 
range [0,1]. We fix the value of the coupling K and we 



let the system to evolve until a stationary (asymptotic) 
value i?oo for the order parameter is obtained. In this 
way it is easy to recognize a Kuramoto-likc transition 
from an incoherent phase (for K < K c ~ 1.4) to a par- 
tially coherent (for K £ [1.4,4.0]) and, finally, to a fully 
synchronized phase (for K > 4.0) 8]. We now focus on 
the details of the dynamical evolution in each of the three 
phases. 

In Fig. 4 we show the case of very small coupling, K = 1.0. 
In the left part we show the time evolution of the opinions 
and of the opinion changing rates (angular velocities or 
frequencies). In the right part, instead, we plot the final 
distribution of opinions and the order parameter time 
evolution. Because of the weak interactions we are in 
the incoherent phase and each agent tends to keep his 
natural opinion changing rate. It follows that the differ- 
ent opinions diverge in time without reaching any kind 
of consensus. In correspondence, the order parameter R 
takes the minimum possible value that, at variance with 
the Kuramoto model, is not zero. We could look at this 
case as to an 'anarchical' society. 

In Fig. 5 we plot the same quantities than before but in 
the case K — 2.0. The coupling is still weak but strong 
enough to give rise to three different clusters of evolv- 
ing opinions, each with a characteristic changing rate: 
the largest number of the agents, representing what we 
could call the " public opinion" , moves with an intermedi- 
ate rate along the opinion axis, but there is a consistent 
group of people remaining behind them and also a group 
of innovative people (quicker in supply new ideas and 
ingenuity). From a political point of view, we could in- 
terpret this situation as a sort of 'bipolarism' with a large 
number of 'centrists'. In this case the order parameter 
is larger than in the previous example, but still less than 
one since the opinion synchronization is only partial. 
Finally, in Fig. 6 we report the case K = 4.0. Here the 
coupling is so strong that all the opinions change at the 
same rate and we observe a single final cluster in the 
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opinion distribution. In this 'dictatorial' society all the 
agents think in the same way and follow the same trends 
and fashions. Although the natural frequencies of the 
agents are - as in the previous examples - different from 
each others, their opinion changing rates rapidly synchro- 
nize (frequency locking) and thus the order parameter R 
reaches a saturation value equal to one. 
Summarizing, it has been found 8] that in order to ensure 
a 'bipolarism' - i.e. an equilibrium between conservative 
and progressist components - a changing society needs 
a level of coupling K strictly included in a narrow win- 
dow (1.5 < K < 2.5) inside the partially synchronized 
phase. Otherwise such an equilibrium will be broken and 
the final result will be anarchy or dictatorial regime. But 
it is worth to observe that these conclusions have been 
obtained for systems with fixed coupling K , simulating 
societies with a stable degree of interconnections among 
their members. Thus it is interesting to explore what 
happens if the coupling is let to increase its value dur- 
ing the dynamics, in order to simulate a society in which 
the interconnections between the agents increase in time, 
due for example to the improvement in transport or in 
communications. 

In Fig. 7 we show the results for a system in which the cou- 
pling is uniformly increased from K — 0.1 to K = 10.1. 
The agents' opinions initially spread freely, and then 
rapidly freeze in a large number of non-interacting clus- 
ters with different changing rates and variable sizes. Ac- 
tually, it results that this particular cluster distribution, 
that could be socially interpreted as a multipolarism, 
cannot be obtained in simulations with a constant cou- 
pling. This could suggest that the increase of interactions 
between the members of a society is crucial to stabilize a 



plurality of different non-interacting clusters of opinions 
(different ideologies, political parties, etc.) typical of a 
multipolar democracy. It seems to suggest also that a 
stable bipolarism is possible only in societies with a fixed 
degree of internal interconnections. 



III. CONCLUSIONS 

In this paper we have shown that even simple opin- 
ion formation models are able to capture many general 
features of real social systems. In the first part we have 
discussed an extension of the scalar opinion dynamics of 
Hegselmann-Krause model to the case in which the opin- 
ion is not just a scalar but a 2-vector with real-valued 
components. We investigated a community where ev- 
erybody talks to everybody by means of Monte Carlo 
simulations of discrete opinion distributions (for squared 
or circular opinion space) and by solving numerically a 
rate equation for a continuum opinion distribution. By 
studying the consensus thresholds we found that the con- 
tinuum case can be considered as the limit of the discrete 
one for a great number of agents. In the second part 
of the paper we have discussed a social variation of the 
Kuramoto model, the so called Opinion Changing Rate 
model (OCR). The concept of 'opinion changing rate' 
transforms the usual approach to opinion consensus into 
a synchronization problem. As for the Kuramoto model, 
the OCR model exhibits a phase transition from an in- 
coherent phase to a synchronized one and shows many 
interesting features with a clear social meaning. 
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